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Jiang et aJ Reply: In the comment by Konig, Gefen, 
and Silva[]J, their main point is to question the ap- 
proximation we used in solving the Green's function. 
They believe that this approximative method is inappro- 
priate to describe spin-flip-related dephasing processes 
caused by intradot interaction, so they believe that our 
conclusions are ill-founded because of that. We agree 
that we indeed used an approximation in calculating the 
Green's function. However, we believe our approxima- 
tion to be reasonable and the main conclusion, namely, 
an intradot electron-electron interaction can not induce 
dephasing, should hold. 2] In particular, we emphasize 
that we have taken the higher-order terms in solving 
the Green's function than appeared in their previous 
publications. 0, 0_We will compare the approximation 
used their papers (J Q with that in our work.Q 

Before we make detailed comparisons, we would first 
make several remarks: 

(i) They and we used the same formulaes to calculate 
the current. Those formulaes were advanced by Meir and 
Wingreen.|| 

(ii) In their device as well as ours, the interaction only 
exists in the quantum dot (QD) and the other parts of 
the device should be non-interacting, so all Green func- 
tions (except G dd (uj)) can be exactly expressed using the 
intradot Green function G dd (uj). This means if G r dd {uj) 
is obtained, then all other Green functions as well the 
current can be calculated straightforwardly without any 
further approximations. 

Here we give an example to exactly express G r sk d (oj) 
and G fk,d( u ) b y G dd {uj) in our system. G r sk d {u) and 
Gij(w) are the Fourier transforms of G r sk d (t) and 
G< M (i) with GJ M (t) = -10(f) < {c sfe (i),4(0)} > and 
Gf k M) 



i < d d (0)c s k(t) >. Following the process of 



Phys. Rev. B 50, 5528 by Jauho et. al,p we have 



Glfc.dM - 9 r s ktsiG r i d + g r sk t s2 G 



2d 



ild 



9 r s k,i^G dd + 
+9 r s kAt±G dd 



9fk,i^G dd 

+ QskA^Gd, 



dd 



where g r,< are the Green functions of the device decou- 
pled to the QD (i.e. with t% = t± = 0) and they can be 
solved exactly. Similarly, all other Green functions can 
also be expressed using G dd (cu) and G dd (uj). Moreover 
although the Keldysh Green function G dd (w) in general 
can not be expressed by G dd (iv), J dujG dd (ui), that is ac- 
tually needed in the calculating current, can be expressed 
by G dd (oj). For example, in their paper, Q they get this 
relation by using II + Jr = [see their Eq.(3.8) and 
(3.11)]. In our paper, Q we can obtain the corresponding 
relation at large T case by using the steady state condi- 
tion. 



(iii) By using those exact relations among the Green 
functions, the current can be expressed solely by the in- 
tradot Green function G 



lj). In their paper, [4j they 
give those expressions [see their Eq.(3.9) and (3.12)] as: 



dd 



AO) _ 

1 u — 



r(i) 



4e TlT} 



h T, 



r, 



duImCrJ®(f L -f R ) 



^V^lTfllWl C0Sl P J dujRe G dd \fL - L 



where I R and I R are the zeroth-order (i° e /) term and 
the first-order flux-dependent (£„/) t erm , respectively. 



G dd is the zeroth-order term of G 7 dd , i.e. G r dd 



T dd 



G^Uef + G r dd 2) tl ef + .... While tref is very small, G d 



dd 

is almost same with G^^ . In our work, we use different 
approach. We calculate G dd , then other Green functions, 
and at last the current. We emphasize that those two 
approaches are essentially the same. 

Now we compare the approximations used in their 
work and ours in calculating the Green function G 7 dd \ 
or G T dd for a very small t re j. 

In our work,[2| we first exactly calculate the isolated 
QD Green functions g dd (u) = [u — tdo — U + Un 9 ]/[(uj — 
ed<r)(w — £dcr — U)]. As a second step, we use the Dyson 
equation to obtain G r dd : G r dd = g r dd + g^hg^hG^ + 

^*4544**^ + ^tl5M*4^ + ^*455l*lG5«I' We a g ree 

that the second step is not exact, but it is a fairly good 
approximation while the QD is weakly coupled to other 
parts of the system. 

Furthermore, if only for analytical results (i) and (ii) 
in the first paragraph of the right column on page 3, not 
for the numerical calculations, we may loosen the above 
approximation. We may first consider that the Green 
functions G r for the system decoupled to the source and 
drain leads has been exactly solved. Secondly, using the 
Dyson equation to obtain the Green function of the whole 
system. Then the results (i) and (ii) on page 3 can also 
be obtained in a straightforward manner although G r is 
still unknown. 

Next, let us examine what approximations are used in 
their papers. 0, Q (i) For U = 0, they get G r } d ] \u) = 
l/(u> — ed + i0 + ), where is the intradot level, (ii) 
For U = oo, in calculating they take G^°' ) = 
(Po + P a )/(w -e d + ,0+) = 1+ / (£d)tJ _J + .o + - (hi) For 
U = oo, in calculating [i.e. to obtain Eq.(3.16) 
from Eq.(3.9) in Ref.0 ] they take G r } d ] = (P Q + P T + 



Pi) 



u>-e d +i{V L +T R )/2 ui-e d +i{T L +T R )/2 ' 



So they do not 



calculate G^ ** at all and they directly write down G dd °^ 
from their intuitive picture. In particular, for U = oo 
they use different expressions of G r dd ^ in the currents 
J(°) and/W. This is a serious error because there is only 

r(0) 

one G dd and it can not be given two different values. 

At last, we reply their other three comments. (1) They 
comment that tt > 1 near resonance in our Fig. 2 inval- 



2 



idates it as a good measure of coherence. In fact, this 
has been emphasized in our Letter, e.g. see the para- 
graph after Eq.(l), or the left column of page 3, etc. We 
mention it here again: If only the first-order tunneling 
process exists, tt describes the degree of coherence; oth- 
erwise when the higher-order tunneling processes are not 
negligible , tt as well as the amplitude of conductance 
Gi does not reflect the degree of coherence. In our Letter 
we design a system (i.e. open multi-terminal AB setup) 
in which the first-order tunneling process dominates, and 
we carry out a study of tt in such a system. (2) They 
comment that our Eq.(4) is wrong, as it relies on the 
single-particle formalism. Notice in Eq.(4) we discuss 
the case of U = [see the paragraph before Eq.(4)]. (3). 
They comment that AG(4>) should be zero at 4> = 0. We 
made an print error in figure capture, AG should be de- 
fined as AG = G(4>) — Go- Here we also show the curves 
(see Fig.l in this reply) for AG = G(<f>) - G{<p = 0). We 
gratefully acknowledge them for pointing this out. 

In conclusion, all our results should hold. The e-e in- 
teraction does not induce any dephasing effect and the 
asymmetric amplitude does not associate with the de- 
phasing effect. 
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FIG. 1: AG = G{4>) - G{4> = 0) vs 4> with the same param- 
eters as the Fig. 2b in Ref.0 ] 
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